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The earliest results of CMS exhibit central pseudo rapidity densities larger than the predictions of the different 
models. Introducing on this basis new guidelines with larger multiplicities of secondaries in the models 
implemented in the simulations, we examine the consequences in 7 ray families (spikes in rapidity distribution, 
coplanar emission) and very large EAS (penetration power in the atmosphere). 
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1. INTRODUCTION 

The collection of high energy cosmic ray events per- 
formed with X ray emulsion chambers (XREC) in the 
energy range covered by the LHC has been the origin 
of questions for high energy physics. Conversely, the 
experience of the colliders and the earliest results of 
the LHC help considerably the interpretation of those 
events by sophisticated simulations; the separation of 
genuine hints of new physics among complex effects of 
natural fluctuations becomes more reliable. Usually, 
the gamma ray families recorded with XREC at very 
high altitudes are not contained events; the origin of 
the collision can be only approached by the invariant 
mass of the 7r 's and sometimes by the geometry of 
the most energetic subcascades. 

Aside from those events which are more or less de- 
pendent on the primary composition, the extensive air 
shower data, through the characteristic properties of 
the different components is also correlated with the 
important features of the primary collisions up to en- 
ergies larger by 3 orders of magnitude (in laboratory 
system) than the present limit of the LHC. Additive 
uncertainties come from the cosmic ray particles in- 
teraction with the air nuclei (N, O, Ar) instead of 
nucleons, the lower energy limits of the experiments 
in p-A collisions and a lack of observation in the frag- 
mentation region in colliders. 




Figure 1: 7's pseudo rapidity distribution in Laboratory 
System for 149 7's recorded in the event JFlafl at 
CONCORDE flight altitude of 17km ; typical spikes may 
be observed near 7.6, 10.3 and 11.4 units of rapidity. 



2. REMARKABLE COSMIC RAY EVENTS 

The Centauros have been widely studied in the 
previous ISVHECRI symposia and measurements are 
done in colliders in relation with iV 7 — N c h relation. 
Hence, this section is concentrated on other striking 
features such as the abnormal fluctuations in pseudo 
rapidity distribution and the coplanar emission. The 
observation of the pseudo rapidity distribution near an 
energy above 10 6 GeV is rare ; we may quote the Texas 
Lone Star, a dozen of A-A collisions recorded in bal- 
loons by the JACEE and an event recorded in the Con- 
corde experiment [l[ exhibiting spikes in the pseudo 
rapidity distribution (fig. 1) . The localisation of the 



vertex in the cabin wall of the Concorde airliner by 
both triangulation and invariant mass method allowed 
the reconstruction of this distribution, also plotted in 
the Centre of Mass System (CMS) on fig. 2 assuming a 
primary energy of 10 6 GeV. The forward cone (142 7's 
over 149) is compared here with the charged inclusive 
and semi inclusive NSD pseudo rapidity distributions 
for Z kno = 2 and Z kno = 3 (Z kno = N ch / < N ch > 
is the Koba-Nielsen-Olesen parameter) averaged over 
3000 events generated with our Monte Carlo generator 
(see section 3). 

The question occurred how those spikes could arise 
of simple random fluctuations or might indicate the 
formation of hot spots of quark-gluon plasma @ , the 
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Figure 2: 7's pseudo rapidity distribution in CMS for 142 
7's (solid line) recorded in the event JFlafl ; the spikes 
are shifted near 2.8 and 3.9 units of rapidity. Simulated 
NSD distributions are superimposed for the inclusive 
component (dotted-dashed)and for semi inclusive data 
with Zkno = 2 (dotted) and Z\- no = 3 (dashed) 



peaks of fig. 1 being considered as the drops of hadronic 
matter into which the quark-gluon-plasma condenses. 

Such events have been also searched in UA5 CERN 
experiment Q : the random grouping of clusters in 
Monte Carlo simulated events suggest that there could 
be no need to include here new phenomena. In the 
case of JFlafl which is a p-Al collision, a similar in- 
vestigation might progress by the analysis of individ- 
ual 7's distributions between Z^ no = 2 and Z^no = 3 
for events corresponding to fig. 2. 

At energies around 10 PeV, alignments of 7's and 
hadrons have been observed, suggesting a coplanar 
emission of the most energetic secondaries. One pair 
of events in addition to a dozen clear events recorded 
at mountain altitude have been collected under a min- 
imal cascading in the low stratosphere, JF2af2, the 
most energetic event during the exposures on the 
Concorde and Strana during Siberian balloon flights. 
JF2af2 is characterised by a near perfect alignment 
of all the cascades with energy exceeding 50 TeV 
(fig. 3). After extensive simulations, we have shown 
that alignments may appear in emulsion chamber and 
result from rare fluctuations of standard collisions [4j . 
However, both events JF2af2 and STRANA with a 
primary energy near ■s/ls) = 5 TeV exhibit a simi- 
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Figure 3: Most energetic secondaries among the 211 7's 
recorded in the event JF2af2. 50% of the visible energy 
^2 E~f (1586 TeV) is contained in the alignment. 



lar visible energy, large multiplicities and require very 
large transverse momenta (around 10 GeV/c). 

The very modest exposure of the chambers (500 
hours with an area of 0.2 m 2 combined with the impor- 
tant reduction of the primary flux after the knee) ob- 
jects the hypothesis on standard fluctuations. There- 
fore, a mechanism involving the valence quarks and 
diquarks, based on relativistic string fragmentation 
has been proposed one pair q — q or qq — qq is 
created when the distance L separating both valence 
quarks exceeds a threshold value. The string fragmen- 
tation tension corresponds to a tension k = (X/2)ira' 
of about 1 GeV/fm, a' being the Regge slope. The 
transverse momentum of the quarks emitted is related 
to the tension by the relation: 



(1) 



The fragmentation of the new strings stretched be- 
tween valence quarks and diquarks and also between 
the partners of the valence diquarks (assumed to be 
strongly bounded) is the origin of the large pis : a 
strong tension of those last strings (n multiplied by 
10 — 15) provides large pis after a violent diquark 
breaking as in fig. 4 for example where < pt > = 11 
GeV/c>, able to reproduce JF2af2. We have here ap- 
plied the Schwinger mechanism 

The fragmentation of such strings would occur 
above an energy threshold of 200 GeV in CMS per 
valence quark. Near this threshold the most energetic 
secondaries are aligned (minimal energy required for 
maximal tension of the string when 3 valence quarks 
are aligned); this could explain why the coplanar emis- 
sion is more likely observed close to 10 7 GeV. Such 
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Figure 4: pt distribution of valence quarks assuming 
valence diquark breaking and Schwinger mechanism. 

configuration is no more necessary at energies exceed- 
ing this threshold, the excitation energy being widely 
sufficient to generate the valence diquark breaking. 

There are some super families at energies above 
the energy of the LHC (Andromeda, Tadjikistan, Fi- 
anit...). Those events present halos, multiple cluster- 
ing, large multiplicities and large p tl but up to now 
no coplanar emission. 
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Figure 5: Charged pseudo rapidity distributions for the 
inclusive NSD component averaged over 3000 collisions 



dau hydrodynamical model (having overestimated the 
conversion in an s - 25 multiplicity of the major part 
of the CM'S available energy) and the prediction of 
Feynman's scaling in favour of s 013 . 



3. SIMULATION AND EAS DATA 
3.1. New Model Guidelines 

Taking the opportunity of the earliest results of the 
LHC given by CMS [6|, we begin to update and im- 
prove the collision generator HDPM implemented in 
CORSIKA and report here the preliminary results ob- 
tained with the non-diffractive code. The parame- 
ters of the main guidelines such as the average central 
pseudo rapidity density < dN c h/drj >, the average 
multiplicity < N c h > have been tuned on the experi- 
mental data as follows for y/{s) > 0.9 TeV: 

< dN ch /dr) >= 0.595s 013 < N ch >= 2.257s 0195 

(2) 

instead of the previous adjustment to UA5 data: 
< dN ch /dr) >= 0.74s 0105 < N ch >= 7.2s 0127 - 7 

Those asympotic tendencies turn to be an interme- 
diate between the behaviour predicted by the Lan- 



Our collision Monte Carlo generator has provided 
in those conditions the pseudo rapidity distributions 
dN/dr) for charged and neutral particles, as shown for 
charged hadrons on fig. 5. The central pseudorapidi- 
tyies are in agreement with CMS, ATLAS and ALICE 
measurements. Other characteristics of the p-p colli- 
sions are shown in Table |TJ The inclusive and semi- 
inclusive data are compared in Fig. 2 at v/(s) = 1.37 
TeV (comparison with JFlaFl). The discrepancy 
with earliest extrapolations based on UA5 data are 
not visible at this energy, but according to Table Q] 
we can expect clear asymptotic tendencies such as 
larger multiplicities, larger violation of KNO scaling, 
larger central rapidity densities and a larger energy 
flow transferred to secondaries (as suggested by the 
inelasticity Ktot in Table HJ). The microscopic mod- 
els with intrinsic parameters tuned to involve those 
effects will produce for proton initiated showers more 
rations and a maximum depth of the e.m. component 
at higher altitude. Consequently, the present values 
proposed for the primary mass are overestimated in 
the knee region and beyond. 
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Table I Characteristics of NSD inclusive data 



Vis) 


0.9TeV 


2.6TeV 


7.0TeV 


< dN ch /drj > 


3.48 


4.47 


5.95 


< N ch > 


33.6 


49.27 


73.86 


< K to t> 


0.59 


0.64 


0.65 


<Pt> 


0.46 


0.49 


0.54 



3.2. Extensive air showers at UHE 

At ultra high energy, the maximum depth X max is 
also an indicator of the characteristics of the primary 
interaction in reason of the relation of the individual 
maximum of the main subcascade with two factors ; 
the vertex of the earliest collisions dependent on the 
inelastic cross section on one hand, the individual en- 
ergy given to the most energetic 7's according to the 
profile of the fragmentation region on the other hand. 
The situation is compared on fig. 5 with the different 
models [7[ together with the data collected by HIRES 
and AUGER. In contrast with HIRES, above an en- 
ergy of 2 — 3 EeV, AUGER suggests that the pene- 
tration power of cosmic particle in atmosphere levels 
off with an asymptotic tendency pointing to the alone 
event and most energetic event of the Fly's Eye. Is 
here a hint of a phase transition to QGP (an energy 
density larger than 30 GeV/Fm 3 can be expected for 
100 EeV collisions)? Is it the signature of a transition 
to a heavy primary composition? We observe that the 
last bins of energy in fig. 6 for AUGER contains hardly 
20 primaries and the hybrid showers selected by the 
surface detector might not reflect the actual situation 
of the primary flux in those bins. 
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Figure 6: Penetration power of the UHE cosmic ray 
particles in the atmosphere. The solid line with a 
question mark corresponds to the tendency suggested by 
AUGER data 



4. CONCLUSION 

The preliminary data coming from the LHC implies 
at least an incremental revision of the models used 
in EAS simulation. The interpretation of the simula- 
tions will turn consequently to lower the primary mass 
composition at the knee level and above. The modest 
changes expected up to 10 s GeV could be more im- 
portant at AUGER level and would favour the inter- 
pretation of HIRES (GZK behaviour and pure proton 
component). 

In the LHC energy range, the conjecture is similar 
to the situation preceding the discovery of the charm, 
proposed in cosmic rays, rejected in a first accelerator 
experiment and finally discovered in further measure- 
ments. Furthermore LHC data will help to clarify 
the role of complex fluctuations in XREC data: up- 
dated XREC experiments (for instance with emulsion 
bricks) could also be performed in parallel with the 
LHC in order to determine the behaviour of the most 
energetic secondaries in fragmentation region. 



References 



[1] J.N. Capdevielle et al, Proc.l9th ICRC , 6, 127, 

Kyoto (1979), Nuovo Cimento 69, 295 (1982), 

Nucl. Phys. B, 75A, 12 (1999) 
[2] A. Faessler,Proc. Nucleus-Nucleus collisions Int. 

Conf., (Michigan, 1982), Nuclear Physics A400, 

578c (1983) 

[3] UA5 A general study of proton-antiproton physics 
at y/la) = 546GeV, G. J. Alner et al, Physics 
Reports, 154, 247 (1987) 

[4] J.N. Capdevielle, R. Attallah and M.C. Talai , J. 
Phys. G, 31, 373-388 (2005) ; Nucl.Phys. B 175- 
176, 137 (2008); Nucl. Phys.B, 196, 395 (2009) 

[5] C.Y. Wong, Introduction to High Energy Heavy- 
Ion Collisions, World Scientific 31, 65 (1994) 

[6] CMS Collaboration, larXiv:1002.062lV 2 [hep-ex] 
(2010) 

[7] J.N. Capdevielle, F. Cohen, B. Szabelska and J. 
Szabelski, J. Phys. G, 36, 075203 (2009) 



C12 



